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On the basis of the sequence of marginal observables the evolution equations of the 
microscopic phase density and its generalizations is discussed. We introduced dual 
BBGKY hierarchy for these microscopic observables and their average values. 
In the space of integrable functions, for initial-value problem of generalized hydro- 
dynamical equation the existence and uniqueness theorem is proved. 
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1 Introduction 

The aim of the present paper is to derive the rigorous solution of the of generalized 
hydro dynamical equation for average value of marginal additive-type observable, known 
as marginal microscopic phase density. We proved the existence and uniqueness theorem 
for the such solution. 

The paper is organized in the following order. 

In section 2 we introduce the definitions used for the description of the average values 
of observables and equivalent approach for the average values of marginal observables. In 
section 3 we give the definition of the microscopic phase density and its generalization. We 
consider the initial value problem for the dual BBGKY hierarchy for marginal observables 
and rewrite this hierarchy in terms of "macroscopic variables." Using the rule for average 
value in the grand canonical ensemble we derive the equations for average values of k-axy 
type of marginal microscopic phase density. In section 4 we consider the nonperturbative 
solution of initial value problem of such equations. This solution is presented over the 
particle clusters whose evolution is governed by the corresponding order cumulants of 
relevant evolution operators. In section 5 we establish that in fact, if initial data is com- 
pletely defined by a average value of marginal additive-type observable, then all possible 
average value of marginal k-type observable of infinitely many particle system at arbitrary 
moment of time can be described within the framework of an average value of marginal 
additive-type observable by the hydrodynamical equation without any approximations. 
The existence and uniqueness theorem of a solution is presented in section 6. 

2 Preliminary facts 

We consider the system of a non-fixed (i.e. arbitrary but finite) number of identical parti- 
cles with unit mass m = 1 in the space M 3 (in the terminology of statistical mechanics it is 
known as nonequilibrium grand canonical ensemble [IS]). Every particle is characterized 
by the phase space coordinates Xi = (qi,Pi), i.e. by the position on the space G M 3 and 
momentum pi G K 3 . 

An observable of this system can be described by the infinite sequence 

A = (A 0l A 1 (xi), ■ ■ ., A n (x 1 , . . .,x n ), . . .) 

on n particles. The average values of observables (mean values or expectation values of 
observables) can written by the following formula 

oo „ 

(A)(t) = (l,D(0)y\A(t),D(0)) = (l^(O))- 1 ^- / d Xl ...dx n A n (t)D n (0) (1) 

n=0 U ' J 

where (1,D(0)) = X^Lo^T / dx% . . . dx n D n (0) is a normalizing factor (grand canonical 
partition function) . The sequence D(0) = (1, Di(0, x{), . . . , D n (0, xi, . . . , x n ), . . .) of prob- 
ability densities of the distribution functions D n (0) at the initial moment of time and the 
sequence of observables A(t) = (A , A\(t, xi), . . . , A n (t, xi, . . . , x n ), . . .) is a solution of 
initial-value problem of the Liouville equation for observables. If A(0) is the sequence of 
continuous functions and D(0) is the sequence of integrable functions, then functional 
exists. 
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An equivalent approach of the description of evolution of many-particle systems, that 
enables to describe systems in the thermodynamic limit, is given by the sequences of 
s-particle (marginal) observables G(t) = [Go, Gi(t, xi), . . . , G s (t, X 1 « . . . . X £ ),...). The 
sequence G(t) is a solution of the initial-value problem of the dual Bogolyubov chain of 
equations (dual BBGKY hierarchy). In that case the average values of observables in 
moment of time t G K are determined by the functional [T5] 

00 1 f 

(A)(t) = (G(t),F(0)) = J2 -7 / dx x . . . dx s G s (t)F s (0), (2) 

where the sequence of marginal observables G(t) in terms of the sequence A(t) are defined 
by the formula 

s (—1)™ s 

G s {t,X!,...,x s ) = ^ j— A s~n(t,Y\{x jl ,...,x jn }), (3) 

where Y = (x\, . . . , x s ), s > 1, and the sequence ^(0) of marginal distribution func- 
tions is defined in terms of the sequence D(0) as follows (nonequilibrium grand canonical 
ensemble) 

00 1 f 

F s (0,x 1 ,...,x s ) = (l,D(0))- 1 J2- / dx s+l ...dx s+n D s+n (0). (4) 

n=0 U ' ^ 

3 Evolution of marginal microscopic phase densities 

Let N(t) = (N® (*),..., (*),...), where (t) = (0, . . . , 0, ivf } (t) , . . . , ivf } (t) , . . . ) , 
k > 1, be the sequence of microscopic phase densities of k-ary type 

n k 

Ni k \t) =NW(t,Z 1 ,...,£ k ;x 1 ,...,x n )= ]T 5(6 - X H (t, x u . . . , x n )), (5) 

where 5 is the Dirac 5-function, £i,...,£fc are the macroscopic variables & = (i^r*) G 
IR 3 x IR 3 . Functions x±, . . . ,x n )}._ v n > k > 1, are the solution of the Cauchy 

problem of the Hamilton equations for n particles with the Hamiltonian 

n 2 n 

/7 »-E|- E *(*-<&). (6) 

i=l i<j=l 

where $(g^ — <£,) is the two interaction potential, and with the initial data x\, . . . , x n . 

Thus, according to fl3]) at initial moment of time t — the sequence of marginal 
observables of fc-ary type microscopic phase densities (J5J) is given as follows 

k k 

G«(0) = (0,...,0, Y[5^-x H ),0,...). (7) 

We introduce the sequence of marg inal observables [TF] G(t) = (G (1) (£), • • • , G (fc) (*)> • • • ) 
of A;-ary type G^(t) = (0, . . . , 0, Gf \t), . . . , G { s k \t), . . . ) . 
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The marginal microscopic phase densities Gi k \t) = Gs(t, £1, . . . , xx, . . . , x s ) of ev- 
ery k-ary type are governed by the initial- value problem of the dual BBGKY hierarchy [18] 

|G?'w=(i:<*.^>-i:<|*(»-*).|j>w)- 

Gi k \t)\ t=0 =Gi k \0), s>k>l. (9) 
j 

Here Y = {x u ...,x s ) and (xx, . . . , V, . . . , x s ) = (sci, . . . , Xj^x, x j+1 , . . . ,x s ) — Y\ Xj. 

Taking into account definition ([7]) initial value problem OH])- ([9]) can be rewritten in 
terms of variables £. Indeed, in terms of variables £i, • • • , Cfe f° r t ne marginal microscopic 
phase densities of k- ary type G^(t) = (0, . . . , 0, G^\t), . . . , G { s k \t), . . . ) we derive 

|^W = (-E<«..^>+E<|*(n-r j ),^>)G?'W + 

i=i »^i=i 



with the initial data 



+ E / ^(^(n-r w ),^Gf +1) W 



(10) 



with the initial data 

s k 

G^\t)\ t=0 = Yi n*(6-^)^, (ii) 

where 1 < r < s and for k = s, the marginal microscopic phase density Gi s \t) is governed 
by the Liouville equation. 

It should be noted that in terms of variables • dual BBGKY hierarchy (1101) 
is represented as the Bogolyubov set of equations with respect to the arity index k > 1, 
while evolution equations (ITU]) have a structure of a sequence of equations with respect 
to the index of number of particles s > k. 

Averaging initial value problem f[TU]) -( TTT]) the evolution equation of average value ([2]) 
for k-ary type microscopic phase density, according to (j2J) has the form [17] 

= " E ("'■ + E (fair, - r,), 

1=1 i+3=l 

+ E/ ^ + i(^-$(r,-r fe+1 ),^-)(G( fc+1 ))(t), (12) 

1=1 

with the initial data 

(G«)(t,ei,...,6)l*=o = (G (fc) )(0), k>l. (13) 
We remark that due to functional (j2J) initial data f)13p are given as the following functions 

(G«)(o,6,...,6) = ^(o,6,---,a), 

where F s (0, £i, . . . , £ s ) is the value of the solution of BBGKY hierarchy at initial moment 
of time at a point £i, . . . , £ s - 
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4 On solution of initial- value problem for the average 
value of marginal microscopic phase density 

To derive the hydrodynamic equation for one-particle function it is necessary to construct 
a non-perturbative solution of Cauchy problem (I12I) - (I13I) in the following form. 

Consider a solution of Cauchy problem (jl2p - (ll3|) that is defined by the expansion 
over the arity index of the microscopic phase density whose evolution are governed by 
the corresponding-order cumulant (semi- invariant) of the evolution operators (fT6l) . con- 
structed in [17] 

• • • ,&) = (U(t)(G)(0)) k (Ci, ■ ■ ■ = 

oo 

= E- / dH k+1 ...dZ k+n % 1+n (-t)(G^)(0), (14) 

n=0 ' ^ 

where G(0) = (0, G^(0), . . . , (7^(0), . . .) is the sequence on integrable functions and if 
n > the operator 

% 1+n (-t) = 2li+„(-t, {Y}, k + l,...,k + n) = (15) 

e H'^iPi-^n^iH 1 ') 

p : {{y},x\y}=u i x l x iC p 

is the (1 + n)£/i-order cumulant of the groups of operators ffTET) . is the sum over 
all possible partitions P of the set {{^}, k + 1, . . . , k + n} into |P| nonempty mutually 
disjoint subsets Xj C {{^}, X \ Y} = {{Y}, k + 1, . . . , k + n}. The set {Y} consists of 
one element of Y = (1, . . . , k), i.e. the set (1, . . . , k) is a connected subset of the partition 
P(|P| = 1). 

If (G^)(0) are integrable functions [19] . series (fT4l converges for small densities. 
In (JTSJ) evolution operators Sk(—i) act as follows 

{Sk(-t)fk) (6, • • • , 6) := A(3i(-t, 6, • • • , 6), • • • , S fc (-t, 6)), (16) 

where functions Hj(£) = (V^(t), Ri(t)) are the solution of the Cauchy problem of the 
Hamilton equations for "macroscopic variables" 

with the initial data 

Ri(0) = n, Vi(0) = Vi, i = l,...,k. 

We remark that using Duhamel formula solution (I14p of Cauchy problem (TT2l)-(IT3l) can 
be represented as the iteration series of BBGKY hierarchy (1121) [17] . 
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5 Hydrodynamical equations for average values of 
observables 

Consider initial data that completely defined by one-particle function (GW)(0,£i), i.e. 

k 

(G)(0) = {I, (G^)(t,^),...,l[(G^)(t,^,...). 

i=l 

Taking into account that (G^)(0, £j) = F 1 (0,^ i ), initial data satisfy the chaos condition. 

For such initial data there is a possibility to reformulate initial value problem f lT2|) - 
f fT3|) as new one for unknown function (G^)(t, £i) with the functional (G^)(t, £i, ■ ■ ■ , | 
(GM)(t)),k>2. 

Functionals (G^ft,^,...,^ | k > 2 are represented by the following 

expansion over the products of one-particle functions (G^)(t) 

(G^}(t,^,...,^\{G^)(t)) : = 

00 1 /" k+n 

= E~ d^ k+l ...d^ +n ^ 1+n {t,{Y},k + l,...,k + n)l[(GW)(t,^), (17) 

n=0 n ' 1=1 

where evolution operators QJi +n (t), n > defined as follows [16] 



9j 1+B (t, {r}, := n \ E(-i) r E • • ■ E 



,, /;-/))] -...-;/!,, — J 



X 



(n - mi - . . . - m r )\ 

r=0 mi=l m r = l v ' 

x2li +n _ mi _..._ mr (t, {Y}, k + l,...,k + n-m 1 -...-m r )x 

r rrij n— mj- rrij+k— 1 fc+n — mi — rfij 

X IIE''' E II (_j _ -J \| x 

j=l r i_g _j _n ij=l V ?i— mj— mj+As+1— »j n—mi—...—mj+k+2—ijJ' 

2 n — m-^—... — m.j-\-k 

x2t 1+ j j (£,« 3 -,fc + n — mi — . . . — ra,- + 1 + 

R— mi — ... — m ^ +fc+l-ij n — m ]_ — . . . — mj + fc + 2 — i ^ 

+ r fc+n-mi-...-m J +2-i j ' ■ ■ ■ , + n — mi — . . . — TOj + nUn-mi-...-mj+1-ij ) , (18) 

where = mj, r J n _ mi _ _ mj+fc+1 = 0. We denoted 2li +n (£) is the (n + l)-th order cumulant 
of scattering operators Sk(t), k > 1 (S\(t) = I is a unit operator) 

i, . . . , k) = s k (-t, i, . . . , k) n Si(t, i). (19) 

i=l 

The samples of two first terms of evolution operator QJ(t) are the following 

9j 1 (t,{r}) = ai(t,{y}), 

k 

9J 2 (f, {y}, A; + 1) = 2t 2 (i, {F}, fc + 1) - {F}) E ^(t, «i, k + 1), 

il=l 
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or in terms of scattering operators 
f0 1 {t,{Y}) = S k {t,Y), 

s 

QJ 2 (f, {Y}, k + 1) = S k+1 (t, Y : k + 1)- S k (t, Y) S 2 (t, j, k + l) + (k — l)S k (t, Y). 

3=1 

On integrable functions the action of scattering operators f)19p are defined by the fol- 
lowing formula 

(S k (t)f k ) • • • , &) = f k (Si (t, E 1 (-t, fx, ... , &)) , . . . , S fc (t, E fc (-t, . . . , £*))), (20) 

where functions B^(t) = (VJ (t) , Ri(t)) ,i = 1, . . . ,n are solutions of the Cauchy problem of 
Hamilton equations for "macroscopic variables" with corresponding initial data. 

The generator of first-order evolution operator ffTgj) . i.e. n = 0, is defined by the 
Poisson bracket with respect to the variables £i, . . . , £ k with an interaction potential $ on 
the continuously different iable functions f k = /fc(£x, • • • )60 

iimi(2j l( t, - i) fk = £ (A $(rj _ ri)j 

In generale case, for n > 1, it holds 

iimi(Qj 1+n (t,{y},x\y))/ fe+n = o. 

For the case of non-interacting particles the evolution operators (fl8|) for n > we have 

limj-(W 1+n (t, {F},X \ r))A+n = <V 

At the same time for initial moment of time t = 

lim ^ (5Ji +n (t, {F}, X \ Y)) f k+n = 5 n , , 

where 5 Uj o is a Kronecker symbol. 

The average value for an additive- type microscopic phase density (G^')(t) is governed 
by the following Cauchy problem (the generalized evolution equation for the average value 
of microscopic phase density) 

|(G«)(t,6) + (^,^)(G'( 1 ))(t,6)= (21) 

= / d&<£*(ri-r 2 ), | {G m m ) 

with the initial data 

(G^)(t^i)\t=o = (G^)(0^ 1 ). (22) 
where functional (G<$) (t,&, 6 | is defined by (HZJ). 
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The functional (G^) (t, 65 6 | (G )(t)) m the collision integral of equation ( 12TT) is 
defined by expansion ( |T7I) 

(G (2) )(t)(t,6,6l(G (1) )(t))= (23) 

= 2ti(t,m)n(G' (1) )(t,6)+ / ^ 3 (2i 2 (t,m,a + i)- 

i=l J 

s 3 

-2Ci(t, {Y}) ^2t 2 (U;,6+i)) HiG^t,^) + 

j=i i=i 

where {y} = (6,6)- 

We represent the first term of expansion (I23"|) in an explicit form. According to the 
definition of two-particle scattering operator (1201) we have 

2 

Mt,{Y})H(G^)(t,&= (24) 



= (GF>)(t, 5i(f, Si(-t, ei,6)))(C (1) )(t, 3a (t, Sa(-i, 6,6))' 
where the functions (i = 1, 2) 

E { (t, ~i(-t, 6,6)) = (V5(-t, 6,6), 6,6) 6,6)) 

are defined as above in (12P1) . 



6 The existence and uniqueness theorem 

Suppose L\ is the Banach space of integrable functions defined on the phase space 1R 3 ™ x 
IR 3 ™, of n-particle system, symmetric under the perturbations of arguments. The norm of 
an element f n of L\ is denoted by 

II = / dxi . . .dx n \f n (xi, . . . ,x n )\, 

L\ C L|, is a subspace of continuously different iable functions with compact supports. 
For functionals of average values of observables (IT7|) the following existence theorem is 
true [T6l 



Theorem 1. //||(G (1) )(0)|| L i (R 3 xR 3) < e~ 10 (l + e- 9 )- 1 then there are exist a global in time 
solution of the initial value problem (fH] ) -( Tf^) of the generalized kinetic equation which is 
defined by 

oo „ n+1 

(G (1) )(t,6) = E- / ^2...rf6 + n2l 1+n (-t)n(C (1) )(0,6), (25) 

n=0 n ' J i=l 

For (GW)(0) G Lj (IR 3 x M 3 ) if zs a strong solution and for arbitrary initial data it is a 
weak solution. 
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Proof. We give a short proof of the theorem here. The full version can be found in |16j . 

Taking into account the fact that operator S(—t) is the isometric operator and for 
\X \ Y\ > 1 the identity 

£ (-1)^- 1 (|P|-1)! = 

P:{X\y}=U l X 1 

is valid, for cumulants 2l(t) we have 

„ k+n 

I d£ k+1 . . . d£ k+n % n (t, k + l,...,k + n) l[(G^}(t, &) = 0, (26) 

Therefore, we can represent marginal functionals of the average value of observables by 
the following renormalized series 



= E~ /^ + l---^ + n^l + n(t,m^+l,---^ + ")n( G(1) )^^)' ( 2? ) 
n=0 n ' ^ 1=1 

where evolution operators 5Ji +n (i), n > defined as follows 



k+n 



2T 1+n (t, {F}, X\Y) := n! £ . . . £ 



;/ /;-/)/ I -...-/llr-l - 



(n - mi - ... - m r )! 

r=0 mi=l m r =l v 1 ry 

x2li + „_ mi _..._ mr (t, {F}, fc + 1, ...,A; + n- mi-...- m r ) x 

x n ' ' ' n / i _ j \\ x 

x2t, , j (t,L,k + n — mi — . . . — rrij + 1 + 

+ri +2 „i . , ■ ■ ■ , k + n - mi - . . . - m 5 + r{ +1 _ { . ) . (28) 
Then we use the following estimate for the cumulants 2t(t) 

» k+n 

J d^...d^ +n \Ql l+n (-t,{Y},k + l,... J k + n)l[(G^}(t,^)\ <n\e n+2 \\(G^)(t)\\ k L t\ 
Using the last estimate for ( 1271) we obtain 

OO 

||(G«)(t|(G«))|| L : < ||(G«)(t)||^e 2 ^||(G'( 1 )>(t)||2, e ri + 

n=0 

oo 

+ 11 (G {1) ) (t) \\ k L ie 3k+2 J2 1 1 (*) I l^ 3fc+2 >\ (29) 



71=1 
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Hence the sequence of marg inal functional of the state • • • , &l k > 2 

exists and represents by converged series ( JT7|) provided that 

|KG«>(t)||^<e-^. 
On the other hand, using the estimate for series ( 1251) . we have 

oo oo 

\\(G (1) )(t)\\ Ll < 1 1 <GC^ > (0) i Li e 2 53 i K G^) > (t) 1 1 2i e 2 ^ e- 1 1 (C^ 1 ) > (0) i |^ - (30) 

n=0 n=0 

Finally, according to estimates ( 129]) and (130]) . we conclude that the collision integral of 
generalized hydrodynamical equation exists under the condition on initial data 



□ 



ll(G (1) )(0)|| L j (R a ><R 3 ) <e- 10 (l + 
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